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ABSTRACT 

The dynamical friction force experienced by a body moving at relativistic speed in a gaseous 
medium is examined. This force, which arises due to the gravitational interaction of the body 
with its own gravitationally-induced wake, is calculated for straight-line and circular motion, 
generalizing previous results by several authors. Possible applications to the study of extreme 
mass-ratio inspirals around strongly-accreting supermassive black holes are suggested. 

Key words: dynamical friction - hydrodynamic drag - extreme mass-ratio inspirals - gravi- 
tational waves 
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The mechanism of dynamical friction (DF), which arises because 
of the gravitational interaction between a massive perturber mov- 
ing in a medium and its own gravit ationally-induced wake , was first 
studied in collisionless systems by Chandrasekhai] i 19431) . and has 
had widespread applications in astrophysics [e.g., stars moving in 
clusters or galaxi es, globular clusters in gala xies, galaxies in galaxy 
clusters, etc.: see lBinnev & Tremainel <1987) section 7.1, and refer- 
ences therein]. In particular, the Newtonian dynamical friction drag 
acting on a perturber of gravitational mass M moving with velocity 
vm in a collisionless system of "particles" with gravitational mass 
m and isotropic velocity dis tribution f(v m ) = dN/(d 3 xd 3 v m ) is 
given ICh^drasekhar 1943) by 



F DF = -16n 2 G 2 M(M + m 



m Io M /(«m)«m^ 
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where A 



,v 2 yp /[G(M + m)], b n 



(1) 

and vtyp being re- 



spectively the maximum impact parameter and the typical ve- 
locity of the partic l es wi th respect to the perturber [see also 
iBinnev & Tremainel dl987l) section 7.1 for a derivation]. The in- 
tuitive reason for the presence of this drag is the fact that the par- 
ticles are attracted by the perturber, which in the meantime moves: 
the particles therefore build up a slight density enhancement be- 
hind it (the wake). It is the gravitational attraction of the wake 
that pulls the perturber back. Note that in the case of a perturber 
moving through a collisionless fluid, dynamical friction is essen- 
tially the only drag force acting on the perturber, besides that due 
to gravitational-wave emission [which is usually s ubdominant be- 
cause it appear s at 2.5 Post-Newt onian order: see IWalker & Willi 
dl980h . and also lPati & WililbOOOl) and references therein] and that 
due to capture of particles by the perturber. 

Equation l[T) was improved by including Post-Newtonian cor- 
rections in iLed ([1969), and was generalized to r elativistic v eloci- 
ties, although only in the weak scattering limit, in lSveil ( fl994l) . The 
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case of a collimated flow of collisionless particles of gravitational 
mass m moving at relativistic speed and impacting o n a perturber 
of gra vitational mass M was instead worked out by IPetrich et al.l 
( 1989) [equation (B17)], who found that in the rest frame of the 
perturber the 3-momentum change is given by 



4itmn oc G 2 M 2 >y 2 [l + (v/c) 2 } 2 



In 



(2) 

where is the number density of particles in the flow far 
away from the perturber and before deflection, v and 7 = [1 — 
(v/c) 2 ] -1 / 2 are the velocity and the Lorentz factor of the perturber 
relative to the flow, and b m i n is the size of the perturber or the cap- 
ture impact parameter 6 m i n ~ 2A/(1 + v 2 )/v 2 if this is a black 
hole. 

Dynamical friction acts also in collisional fluids, together with 
the two other effects mentioned above for collisionless systems 
(i.e., gravitational- wave emission and accretion onto the perturber) 
and ordinary viscous forces, which are not present if the perturber is 
a black hole but instead act if the perturber is a sta£]. However, un- 
like in the collisionless case, it presents different features depend- 
ing on the Mach number of the perturber. The correct behaviour in 
the supersonic case has long been recognized by several authors: 
the steady state Newtonian drag on a perturber M moving on a 
straight-line with velocity v relative to a homogeneous fluid with 
rest-mass density pp and sound spee d c s = v/M (M > 1) was 
found by Rephaeli & Salpetet] d 1980h and bv lRuderman & Spiegel 



1 The drag due to ordinary viscosity is given, for non-relativistic velocities 
and in the laminar regime, by Stokes' law: -Fstokcs = — 6 n V av > a being 
the radius of the perturb er and r) the viscosity co efficient. For instance, in 
a thin accretion disc jShakura"& Sunvaevl[l973l) one has r\ = apoc s H, 
where pp and c s are the rest-mass density and the sound v elocity in the disc, 
H is its height and a ~ 0. 1 - 0.4 tang. Pringle, & Liviol2007l) . Note that 
this drag can be calculated independently of the dynamical friction effects 
considered in this paper. 
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dl97ll) tobe 
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where the maximum impact parameter b max is the Jeans length (or 
the size of the medium if this is smaller than the Jeans length). 
Note that the dynamical friction drag given by equation {3j is com- 
parable to the drag due to Bondi accretion onto the perturber: 
the latter is in fact given b y Fsondi = —vM, where M = 
4\tyG 2 M 2 Po /(v 2 +c 2 ) 3/2 bondi & Hoylel 1 1944 iBondil 1 19521) . 
A being a parameter of order unity. 

Equation ([3} was confirmed by lOstrikd dl999h with a finite- 
time analysis and was generalized to the relativistic case by 
IPetrich et alj dl989l) . who found that equation ((2} remains valid also 
in the collisional supersonic case if the rest-mass density rnn^ is 
replaced by p + p, p and p being the pressure and energy-density of 
the fluid. The physical reason for the presence of a non-zero drag 
in the supersonic case is the fact that sound waves can propagate 
only downwind, inside the Mach cone, producing a non-symmetric 
pattern of density perturbations, which gives rise to a drag by grav- 
itational interaction. 

The subsonic case proved instead to be more elusive. Because 
sound waves can propagate both downwind and upwind, the drag 
is expected to be lower tha n in the supersonic case. In particular, 
Reph aeli & SalpeteJ (l 980) in the Newtonian case and lPetrich et al.l 
dl989h in the relativistic one argued that the drag should be ex- 
actly zero for subsonic motion in a homogeneous fluid, because of 
the upwind-downwind symmetry of the stationary solution for the 
density perturbations excited by the perturber. However, although 
a zero drag can be a useful approximation in many cases, this re- 
sult doesjiot rigorously hold if one performs a finite-time analy- 
sis dOstrikei| [l999). In fact, if the perturber is formed at t = 
and moves at non-relativistic subsonic speed on a straight-line in 
a homogeneous fluid, the d ensity perturbations are give n by the 
stationary solution found by iRephaeli & Salpeten dl980h only in- 
side a sphere of radius c s t centered on the initial position of the 
perturber, and are instead exactly zero (because of causality) out- 
side. The upwind-downwind symmetry of the stationary solution is 
theref ore broke n and the perturber experiences a finite drag, which 
reads (Ostriker 1999) 
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as long as (c s + v)t is smaller than the size of the 
medium. This result was confirmed b y numerical simulations 
( Sanchez-S alcedo & Brande nburg 1999) and was extended to the 
case of a perturber moving at non -relativistic speed o n a circular 
orbit in a homogeneous medium bv lKim & Kin] d2007l) . In particu- 
lar. lKim & KirrJd2007l) found that in the subsonic case the perturber 
experiences a tangential drag, given roughly by equation and a 
drag in the radial direction (towards the center of the orbit), whose 
contribution to the orbital decay is however subdominant with re- 
spect to the tangential drag. Similarly, in the supersonic case the 
tangential drag is roughly given by equation 10 with 6 ma x equal to 
the orbital radius, while a radial drag is present but again remains 
subdominant with respect to the tangential one. 

The purpose of this paper is to g eneralize t o the relativis- 
tic case the finite- time drag found by lOstriken d 19991) and by 
iKim & Kiml d2007l) . While a Newtonian treatment of dynamical 
friction is satisfactory in many astrophysical scenarios, relativis- 
tic expressions are needed in order to study the interaction of 
solar-mass compact objects or black holes with the gaseous mat- 



ter (e.g., an accretion disc) which could be present in the vicin- 
ity of a supermassive black hole (SMBH), where orbital veloc- 
ities close to that of light are reached. Th ese systems, known 
as ext reme mass-ratio inspirals [EMRIs; see lAmaro-Seoane et al.l 
(2007) for a review], are expected to be among the most inter- 
esting sources of gravitational waves for the Laser Interferome- 
ter Space Antenna (LISA), and considerable effort has been spent 
trying to understand whether different kinds of accretion disc, 
when present, can produce an observable signature in the emitted 
gravitational-wave signal. In a series of papers, Karas, Subr and 
Vokrouhlicky c onsidered the interaction between stellar satellites 
and t h in discs dVokrouhlickv ■& Karaslll993l. Il998l ISubr & Karas! 
1 19991 ; iKaras & Subn 1 l200ll) . iKaras & Subn 1 d200ll) , in particular, 
found that the effect of star-disc interaction on EMRIs dominates 
over the effect of the loss of energy and angular momentum through 
gravitational waves in the case of thin discs, both for non-equatorial 
orbits crossing the disc only twice per revolution and for equato- 
rial orbits embedded in the disc, unless the orbiter is very com- 
pact (a neutron star or a black hole) or the disc has a low den- 
sity (e.g., in the region close to the central SMBH if the flow be- 
co mes advection-d ominated). These results agree with those found 
by iNaravanl d2000h . He focused on Advec tion Dominated Accre- 
tions Flows [ADAFs dNaravan & Yil [19941) 1. which were then be- 
lieved to describe accretion onto "normal" galactic nuclejf] (i.e., 
ones much dimmer than Active Galactic Nuclei such as quasars, 
Seyfert galaxies, etc.). He found that for compact objects and white 
dwarfs the effect of the drag exerted by the accreting gas is neg- 
ligible compared with the loss of energy and angular momentum 
through gravitational wa ves, whereas it is not ne gligible for main 
sequence and giant stars. IChakrabartild 199311 19961) studied instead 
the orbital evolution of black hole satellites on circular equato- 
rial orbits embedded in a disc with a non-Keplerian distribution 
of angular momentum, and found that the exchange of angular 
momentum between the disc and the satellite can lead to impor- 
tant effects which have to be taken into account when interpreting 
gravitational-wave signals from such systems. 

Although we do not expect our results to change significantly 
the picture outlined above for ADAFs or ADIOS's, whose den- 
sity is too low to make the effect of the hydrodynamic drag, and 
of dynamical friction in particular, comparable to the effects of 
gravitational-wave emission even if one includes relativistic cor- 
rections, we think that our relativistic corrections could play a more 
important role, under certain circumstances, for black holes or com- 
pact objects moving in higher density environments (Active Galac- 
tic Nuclei, quasars, Seyfert galaxies, etc.). In a subsequent paper 
we will apply ou r results to the case of an accretion flow with a 
toroidal structure (Barausse 2 0071) . 

While our results rigorously apply only to a non self- 
gravitating fluid in either a flat background spacetime (in the case of 
straight-line motion) or the weak field region of a curved spacetime 
(in the case of circular motion), and additional work may be needed 
in order to evaluate the effect of a curved backround, we argue that 
such an approximation is suitable at leas t for a prelimin ary study 
of dynamical friction effects on EMRIs dBaraussell2007l) . Indeed, 
for many purposes a similar approximation is adequate to study 



2 Nowadays, accretion onto "normal" galactic nuclei is believed to be bet- 
ter described by Advection Do minated Inflow Outflow Solutions (ADIOS) 
telandford & BegelmaJ 19991) . However, this is not expected to change sig- 
nificantly Narayan's results because ADIOS's, like ADAFs, are expected to 
present very low densities in the vicinity the central SMBH. 
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gravitational-wave emission by EMRIs around a Kerr SMBH: the 
flat-spacetime quadrupole formula, combined with geodetic motion 
for the solar-mass satellite, gives results which are in surprisingly 
good agreement with rigorously computed waveforms ( Babak et al. 

l2006h . 

Our analysis closely follows that of lOstrikej J 19991) and 
iKim & Kiml d2007l) . and we find that their results still hold for rel- 
ativistic velocities provided that the rest-mass density appearing in 
the Newtonian formulae is replaced by p + p (p and p being the 
pressure and energy-density of the fluid), and a multiplicative fac- 
tor is included: 7 2 [1 + (v/c) 2 ] 2 in the straight-line motion case 
and for the tangential component of the drag in the circular motion 
case; 7 3 [1 + (v/c) 2 ] 2 for the radial component of the drag in the 
circular motion case. 

Throughout the rest of the paper we used units in which G = 
c = 1. 



2 EQUATIONS AND VARIABLES 



fast as r — » oo so as to make the Laplacian operator V 2 invert- 
ible. As an example, let us consider the case of equation (8). First, 
calculating 9 l Tb; and using equation ^ one immediately obtains 



S w = v- a (s i T«; 

Si~ = Toi — 



(10) 
(11) 



where V denotes the inverse of the Laplacian V . Summing over 
the spatial indices trivially gives 



T = <f J Ty , (12) 



and calculating d l d 3 Tij and tFTy using equation ^ one easily 
obtains 



T 



T,t = 2V" 2 d J T l:i - -diT 



(13) 
(14) 



Let us consider a perturber with gravitational mass M, formed 
at t = and moving in a perfect fluid at rest at the ini- 
tial position of the perturbed and having energy density p and 
pressure p there. We write the metric as a Minkowski back- 
ground plus some perturbations produced by the presence of 
the fluid and the pert urber: the gene r al form of such a met- 
ric is known to be teardeerj 1198(1 iKodama & Sasakil 1 1984 
Mukhanov, Feld man. & Brandenbergerlll992l : [Flanagan & Hughes! 
2005b 



ds 2 — g^dx^ 'dx v = — (1 + 2<f>) dt 2 + 2uiidx' dt+ 

fetl-^+^W, X\ = 0, (5) 

where the 3-vector uii can be decomposed into a gradient and a 
divergence-free part, 



LUi = diUJ + LUi 



O LUi 



o, 



(6) 



while the traceless 3-tensor \ij can De split in a gradient part, a 
divergence-free vector part and a (gauge invariant) transverse pure- 
tensor part, 

1 



Xij = DijX + d(iXj) + Xij , D i3 



didj 



di _L at \ 
Xi = o Xij = Xi 



0. 



(7) 



where V 2 = 5 IJ didj. Note that Latin indices are raised and low- 
ered with the Kronecker delta 8ij. Similarly, the stress-energy ten- 
sor can be written as 



T^dx^dx" = Tudt 2 + 2(diS n + S^dtdx* 
T , 



+ 



3^ 



dx^x J , (8) 



where 



0. 



(9) 



Note that the decompositions outlined in equations l[5) and l[8} are 
defined unambiguously if the perturbations go to zero sufficiently 



Inserting equations J 1 2b . J 1 3b and d!4t into equation JS}, one can 
finally derive an explicit expression for the gauge invariant trans- 
verse traceless perturbation E J . Similar considerations apply to the 
decomposition (O of the metric. 

We should mention that our perturbative expansion relies on 
two parameters ei and ei. Deviations of the metric away from a flat 
background are due to the presence of the fluid, which causes per- 
turbations of dimensionless order < ei = 0(C/Xj) 2 [C being the 
characteristic size of the medium and Aj = c s /(47r(p + p)) 
being a generalized Jeans length], and due to the presence of 
the perturber, which is expected to cause perturbations of order 
£2 = M/r, where r is the distance from the perturber. Note that 
the perturbations of the first kind are small if the fluid is not self- 
gravitating (i.e. if C <C A j), while those of the second kind in prin- 
ciple diverge if we consider a point-like perturber. In order to retain 
the validity of the perturbative expansion, we therefore have to in- 
troduce a cutoff r m j n , which is taken to be the size of the star acting 
as the perturber or, in the case where the perturber is instead a black 
hole, the "capture" impact parameter r m in ~ 2M(l + v 2 )/v 2 (i.e., 
the impact parameter for which a test-particle is deflected by an 
angle ~ 1 by the black hole). This ensures that £2 is small and 
can be treated as an expansion parameter. The gravitational field 
produced by the perturber on scales smaller than the cutoff gives 
rise, when coupled to the fluid, to accretion onto the perturber. This 
gives additional contributions to the drag, but these effects can eas- 
ily be calculated separately: see for instance IPetrich et al.l dl989h 
[equation (2.40)] for the drag-force due to accretion onto a black 
hole. When acting directly on the perturber, the gravitational field 
produced by t he perturber it self gives rise instead to the so-called 
self-force [see lPoissonl d2004h for a review], the dissipative part of 
which accounts for the energy and angular momentum lost through 
gravitational waves. 

In order to exploit as much as pos sible the calculations don e in 
the Newtonian case bv lOstrikerj \ 19991) andlKim & Kiml (120071). le t 
us choose the so-called Poisson gauge dMa& Bertschi ngeril 19950 . 
defined by the conditions diUi 1 = diX %1 = 0. In this gauge the 
perturbed metric is 



i Note that this condition can always be satisfied by performing a suitable 
boost. 



ds 2 = - (1 + 20) dt 2 + 2ioidx i dt+ 



\5 ij (l-2i>)+xl 3 ]dx i dx 3 , (15) 
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and the linearized Einstein equations give 



V 2 t/> = 


4vrT tt , 


(16) 


d t ij) = 


4ttS' 11 , 


(17) 


V Ui = 


— 16ivSi > 


(18) 


V 2 = MTtt 


+ T) - 35?^ , 


(19) 




= 87rE" , 


(20) 


dtuJi = 


-87rE+ , 


(21) 


D x5 = - 


-167tEJ , 


(22) 



stein equations one gets the following relations between the matter 
fields: 



V 2 S" 



d t T tt 



V 2 E" = ~(3dtS {l - 
V 2 E, X = 2d t Si 



T). 



(23) 
(24) 
(25) 



which can also be derived directly from the conservation (to first 
order) of the stress-energy tensor with respect to the background 
0. 

Let us now write the stress energy tensor as T M „ — j_ 



metric, d^T^ 



Tifluid . 



T£* Tt . The stress-energy tensor of the fluid is 



rpn 
■*■ it 



(26) 



- — (p + P)u^u v + pg^ , 

where the perturbed metric is given by equation M5\ and p, p 
and fi M are the perturbed energy density, pressure and 4- velocity of 
the fluid: 



P = 



P + Sp . 

= 8m , 



P + Sp, 
-l-<t> 



(27) 
(28) 



(the equation for u, comes about because the fluid is at rest at the 
initial position of the perturber, while the equation for u± follows 
from the normalization condition g^u^u" = — 1). The stre ss- 
energy tensor of the perturber is [see for instancejPoisso^ ^20041) 1 



Tp H (x,t) = M: 



.pert 



(*)), (29) 



''pert i 



where u£ ert 



and x pclt (t) are the perturbed 4-velocity and spatial 
trajectory of the perturber, which for mathematical purposes is con- 
sidered to be a point-particle, and g is the determinant of the per- 
turbed metric U5\ . If one expands the trajectory x pcrt and the 4- 
velocity Sp ert of the perturber as the sum of their unperturbed val- 
ues x pclt and u[J crt plus some perturbations due to the presence 
of the fluid (and therefore of order ei) and some perturbations due 
to the interaction of the perturber with its own gravitational field 
(and therefore of order £2), and notes that g — — (1 + 20 — 6ip) = 
— 1 + 0(ei, £2), equation i29\ can be written as 



T»r(x,t)=M- 



-S (3) (x-; 



.pert 



(«)) 



-'pert 



x[l + 0(si,s 2 )]. (30) 



Note that because of the presence of the factor M = rE2, the stress- 
energy tensor T p ° rt is an intrinsically first order quantity, and drop- 
ping the second order terms, as we have done earlier, we can simply 
write 



T*?\x,t) = M- 



5^(x- 



pert 



(*))• 



(3D 



Perturbing the expression for the conservation of the baryon num- 

(n = n + Sn being the 



ber of the fluid, <9 M ((-g) 1/2 mt M 



perturbed number density), one gets 
8n 



dt 



+ di8u — Sdtip — , 



(32) 



whereas perturbing the Euler equation 5 M = u a V a u M = — (g A " / + 
u' 1 u v )d v h/h [h = (p + p) /n — h + Sh is the perturbed specific 
enthalpy] one obtains 

8n 



d t 5u + di(f> + d t iOi + c s di 



0. 



(33) 



where c s = (dp/dp) 1 ' 2 is the velocity of sound and where we have 
used the first law of thermodynamics (8h/h = c 2 8n/n). Combin- 
ing the divergence of equation d33t and the time derivative of equa- 
tion d32t and finally using equation dl9t , one gets the following 
wave equation for the baryon density perturbations: 

(9? - c 2 V 2 )— = V 2 + 39 t V = 4jt(T« + T) . (34) 
n 

In the next sections we will solve the wave equation J34t and 
the Einstein equations j!6 M22t . It is understood, however, that 
since these equations are linearized, the solutions that we find "in- 
trinsically" have a relative error C(ei, £2)- This error is not to be 
confused with those which we will introduce when solving these 
equations approximately. We will explicitly keep track of the latter 
in the next sections, while we will re-introduce the relative error 
C(£i,£2) due to the linearization procedure only in the final re- 
sults. 



3 STRAIGHT-LINE MOTION 

Let us first consider the case of a perturber moving along a straight- 
line, which is taken to be the z-axis of a Cartesian coordinate 
system: the unperturbed trajectory of the perturber is therefore 
x pclt (t) = y pett (t) = 0, z pclt {t) = vt and the unperturbed 
4-velocity reads Up Crt 9/9a; M = j(d/dt + vd/dz), with 7 2 = 
1/(1 — v' 2 ). Denoting by H(t) the step function, equation d34t can 
be rewritten as 

(d 2 t - c 2 V 2 )— = 4nMy(l + v 2 )5(x)S(y)S(z - vt)H(t) 
n 

+ 4tt[p(1 + 2<f)) + 3p(l - 2V>)] + 4tt(1 + 3c 2 )(p + p)— . 

n 

(35) 

Solving this equation is complicated by the presence of the terms 
4?r[p(l + 2(f>) + 3p(l - 2i/>)] and 4tt(1 + 3c 2 )(p + p)Sn/n on the 
right-hand side. If these terms were not present, we could simply 
solve equation d35 b by using the Green's func tion of the flat w ave 
operator — d 2 + c 2 V 2 , and proceeding as in lOstrikeri (l999) we 
would get 

M-y(l + v 2 ) 



-''pert 



n ( X,t ) * ? c1[(z-vtr + (.,--' -:- '/-)(.!. - M-)}> 1 
where A4 — v/c s is the Mach number and 

1 ifx 2 + y 2 + z 2 < (est) 2 



2 ifM> 1, x 2 + y 2 + z 2 > (c s t) 2 , 
(z - vt)/y/x 2 +y 2 < ~VM 2 - 1 
and z > c s t/M 

otherwise 



(36) 



(37) 
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Note that performing a boost to the reference frame comoving with 
the perturber (the "primed" frame) this "approximate" solution be- 



If u)/c s 2> 1/Aj, one can series expand equation J4U and get 



Sn _ M7 2 (1 + d 2 ) 
n cfr Vl -.M 2 sin 2 6' 



AT 



1 - C s 2 2 



(38) 



where r' and 6 1 ' are the radius and polar angle in the primed frame 
i.e., r' — \J x' 2 + y' 2 + z' 2 and cos 9' — z'/r' in terms of the 
Cartesian coordi nates x', y' and z 1 . E quation J3 8fc agrees with the 
solution found in lPetrich et alj dl989h [equati on (B30)], except for 
the different value of / [this happens because IPetrich et al.l dl989h 
considered the s tationary solutio n instead of performing a finite- 
time analysis: rf lOstrikeJ ( fl999l) for more details]. 

It is not difficult to see that equations d36l l and l |38t are actually 
approximate solutions to equation J35 b - Indeed, the term 47r(p+3p) 
on the right-hand side of equation {35} simply gives rise to an error 
0(C/Xj) 2 in the solution. This error represents the correction due 
to the fact that having a fluid with constant p and p together with the 
Minkowski metric is not a solution of the Einstein equations. When 
it comes to the term 8n (p<j) + 3pip), let us note that the gravitational 
potentials (f> and ip consist of a part of order 0(C/\j) 2 due to 
the presence of the fluid and a part of order — M/r ~ —Sn/n (r 
being the distance from the perturber) due to the presence of the 
perturber. The first part of the potentials therefore gives rise, when 
inserted into the term 8n (p<j> + 3ptp), to an error much smaller than 
the error 0(C/Xj) 2 coming from the term Aiv(p + 3p). The second 
part of the potentials, when inserted into the term 8n(ptp + 3pV)> 
gives rise instead to a Yukawa-like term similar to the term 4n(l + 
3c 2 )(p + p)8n/n appearing on the right-hand side of equation <35t . 
It is not difficult to see that these Yukawa-like terms give rise to a 
relative error ej ~ 0(C/Xj). To see this, one can simply Fourier- 
transform equation l |35l > with respect to time in order to get rid of 
the time derivatives. One is then left with an equation of the form 



5n 



S(x,u) + (a/\ 2 -u; 2 /c 2 ) — (x,u;), (39) 
n 



where for simplicity we have used the same symbol for Sn/n 
and its Fourier transform with respect to time, a is a constant and 
S(x,ui) is a suitably defined source function [inspection of equa- 
tion ( 1351 ) actually reveals that S(x,u>) ~ exp(zo;z/i>)]. Using the 
Green's function of the Yukawa operator V 2 — p 2 (p, being a con- 
stant)Q 



Sn 
n 



(42) 



,3 / f exp [iac s \x - x'\/ (2u)X 3 )] 
A-k\x — x'\ 



Jd 3 x'{'- 



x exp(— iuj/c s \x — x'\)S(x' * 
3 ,exp(-iuj/c s \x - x'\)S(x' ,lu) f c s ej 

x \ — i 71 X I 1 H — 

Atyix — x'\ \ luAj 



and from the last line of this equation it is clear that one gets the 
solution which would have been obtained by neglecting the term 
a/A 2 x Sn/n in equation d39ll, with a relative error c„e.t /(m\i) <C 
£j. For frequencies u)/c s <C 1/Aj (i.e., for wavelengths larger than 
the generalized Jeans length Aj) this procedure is not applicable. 
However, it is clear that for ui — equation J4U becomes the 
solution which would have been obtained by neglecting the term 
a/A 2 x Sn/n in equation l !39t , corrected by a factor ~ (1 + ej). 
Moreover, because S(x, ui) ~ exp(icuz/v), the integral appearing 
in equation J41t averages out if u) ^> v/C. Therefore, the spectrum 
of Sn/n extends up to ttWoff ~ v/C, and the effect of the frequen- 
cies ui/c s <C 1/Aj on the final solution 5n/n(x,t) is negligible 
because Cutoff ~ v/C S> 1/Aj if the fluid is not self-gravitating. 

As such, since we already know the solution of equation d35b 
if we neglect the terms 47r[p(l + 20) + 3p(l — 2ip)\ and 47r(l + 
3c 2 )(p + p)8n/n on the right-hand side [equation J36H. we can 
write the following approximate solution for equation J35t : 

S J1 (X t) = f M^l + v 2 ) 

n V ' ' J c 2 [(z-vt) 2 + (x 2 +y 2 )(l-M 2 )] 1 / 2 

x {l + ej)+0(C/Xj) 2 , (43) 

where, as explained above, the error 0(C/\j) 2 comes from the 
term 4n(p + 3p) on the right-hand side of equation d35t . while 
the error ej comes from the terms 47r(l + 3c 2 )(p + p)8n/n and 
8-k{p4> + 3p4>). [Note that £j(x,y, z,t) — ej(—x,y,z,t) and 
£j{x,y,z,t) = ej(x,—y,z,t) due to the cylindrical symmetry 
of the problem.] Both of these errors are negligible if C <C Aj (i.e., 
if the fluid is not self-gravitating). 

The trajectory of the perturber is governed by the geodesic 
equation of the physical, perturbed spacetime (i.e., the one with 
metric g M „ = ?7 M „ + <5g M „). The familiar form of this equation is 



d 2 x^ 

+ La " df 



df 2 



df a Hf 13 

dT~ 



= o, 



(44) 



G(x) 



exp (— p\x\) 
A-k\x\ 



this equation can be solved and gives 



(40) 



where :Ep Crt and f are the perturbed trajectory and proper time 
while the T's are the Christoffel symbols of the perturbed space- 
time. This equation can be easily expressed in terms of the back- 
ground proper time r, 



Sn 
n 



I 



^pi-^a/Xj-^Mlx-x']) 

4-1TX-X' 



(41) 



4 Using spherical coordinates and the fact that V 2 (l/|a;|) = 
— 47T(5( 3 ' (a;), it is indeed easy to check that (V 2 —fj, 2 )G(x) = —5^(x). 



u ■''pert u " L pert ""^pert 



dr 2 



dr dr 



which can be also written as 



rf 2 r' i dr a dr 13 

11 *^pert . p^t "'■^pert ""''pert 



dr 2 ' dr dr 

dx^Jdr d 

y Q /3 dT dT 



d 2 T (df 

"df 2 Vdr 



dx" 



dr 



(45) 



^pert ^pert ,„,, 

dr dr ■ m 



Using now 



v + Sg^u, equation l l46l can be easily 
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rewritten, to first or der [i.e., neg lecting as usual errors of order 
C(s?,el, £i£ 2 )], as jPoissonll2004l) 



pert 
1 

2 



u x pcrt 

dr 2 



(47) 



The metric perturbations <5g MI / appearing on the right-hand side of 
equation ( |47| > consist of a part produced by the stress-energy of the 
fluid (Sg F v ) and one produced by the perturber (5g F v ). The latter 
contribution, as already mentioned, gives rise to accretion onto the 
perturber and to the self-force. The drag due to accretion is easy 
to calculate separately, as mentioned previously, w hile the self- 
force is in general hard to deal with ( |PoissorJl2004T) . However, it 
is well-known that the self-force is zero in a Minkowski spacetime 
for geodetic (i.e., straight-line) motion in the Lorenz gauge. Since 
the right-hand side of equation J47b is not gauge-invaria nt, the self- 
force itself is not gauge-invariant tearack & Orill200ll). Neverthe- 
less, it is possible to show that at least the dissipative part of the 
self-force (i.e., the one accounting for the deceleration due to the 
loss of energy and angular momentum through gravitational waves) 
is gau ge-invariant and therefore zero also in the gauge which we are 
using dMinoll2003h . (Alternatively, this can be understood from the 
fact that a perturber moving on a straight-line does not emit energy 
through gravitational waves in the quadrupole approximation.) It 
should be noted that the presence of the fluid does not alter these re- 
sults. In fact, one can insert the decomposition Sg^ = Sg F „ + Sg F v 
into the Einstein equations, and split them into equations for Sg F v 
and equations for 5g F v by including in the right-hand sides of the 
equations for t>g F „ only quantities containing the stress-energy of 
the perturber and 8g F v itself. In particular from equations J 1 6b . d!9b 
and d22b . using equation d26b one gets 



■ 47T 



= 8^+4^^, 

P + 2p4> h +(p + p) — 
n 

- J> F = 8ttE! 



V ~~ <V — ""^pert ! 

<p F - <f> F = , 
= -167r(£l pcrt 



n TF 



+ PXij 

F 



(48) 
(49) 

(50) 
(51) 
(52) 
(53) 



From equation {33} it follows instead that Suj 



j^Qand there- 



fore st ama 

with equation d!8b gives 



-(p+p)5u i +pu)i = (2p+p)cjj , which together 



-16ttS 

. F 



X pert 



167r(2p + p)u!i 



167r(2p + pK 



(54) 
(55) 



From equations d48b . d50b , d52b and d54b it therefore follows that the 
metric perturbations Sgf^ produced by the perturber are the same 
as in the absence of the fluid, except for the presence of the terms 
8irp4> F , —l6TvpxTj F and — 167r(2p + p)lo^ f on the right-hand 
sides of equations d48l >. d52b and d54b . Using the Green's function 
of the Yukawa operator it is easy to see that these terms produce a 
contribution of order 0(pr 2 nin d8g F „) ~ O(pM) to the gradients 
d5g F v ~ M/r^i n . To be more specific, let us consider for example 



5 We are making here the simplifying but reasonable assumption that no 
vortical modes Suf- and uif- are excited in the system before the perturber 
is turned on at t = 0. 



the case of ip F . Using equations d40b . d48b and d50b . the solution for 
Tp F reads 



ip (x,t) 



(56) 



3 , exp (-^/8tvp\x - x'\) 



[T%* t (x',t)-WKpT,l elt ( x ',t)] 



Taking now the derivative with respect to x and expanding the ex- 
ponential, it is easy to check that the presence of the fluid simply 
adds a contribution of order 0(pr 2 nin diip F ) to the gradient diip F ■ 
It should be noted that a contribution of order 0(pr 2 nin d5g F „) ~ 
O(pM) to the gradients dSg^^ corresponds to a contribution of 
order 0(p A/ 2 ) to the drag: this contribution can be interpreted, as 
we have mentioned, as being due to accretion onto the perturber. 

We will therefore focus on the force produced by the gravi- 
tational interaction with the fluid, which includes dynamical fric- 
tion. From equations d49t , d5U , J5 3 b and d55t it follows that the 
fluid can only excite the metric perturbations (f> and tjj. Using equa- 
tions d40b . d49b and d51b . we can easily get expressions for the gra- 
dients dn4> F = dfj,ip F evaluated at the position of the perturber 
x — y — 0, z = vt, which enter equation ( 147b . In particular, the 
solution for tp F is 



3 , exp {-y/Mp\x - x'\) 



Sn , 
p + (p + p) — (x ,t) 
n 



(57) 



and taking the derivative with respect to x, one easily gets 



di^ F (x) 



,3 l X - X 
a x 1 r 



P+(p + p) — (x ,t) 



X (1 + EYuk 



(58) 



where we have introduced the error EYukawa ~ 0{pC 2 ) which 
arises when expanding the Yukawa exponential. In particular, note 
that the source p appearing in the integral of equation d58b simply 
gives the gravitational force exerted by the unperturbed medium on 
the perturber. This force is exactly zero if the medium is spherically 
symmetric with respect to the perturber, but in general the net ef- 
fect on the gradients di(f) F = diip F can be non-zero and at most 
of order p C, depending on the shape of the fluid configuration and 
on the position of the perturber. Similarly, the term (p + p)Sn/n 
appearing in equation d58b can be considered as the sum of two 
parts, one coming from the error 0(£/\j) 2 appearing in equation 
d43b and the other one from the rest of this equation. Note that the 
first part is present even if the mass of the perturber goes to zero 
and represents the force exerted by the density perturbations which 
appear because, as mentioned earlier, a fluid with constant p and 
p together with the Minkowski metric is not a solution of the Ein- 
stein equations. The contribution to the gradients <9 M F = 9 /J -!/) F 
from this term can be as large as pC(C/Xj) 2 , and in what follows 
we will group it together with the contribution from the term p ap- 
pearing in equation d58b into a correction e n otDF < ®{p£)- The 
rest of the term (p + p)8n/n gives instead the force exerted by 
the density perturbations produced by the perturber i.e., dynami- 
cal friction. In particular, using equation d43b in equation d58b one 
obtains, for the x and y gradients evaluated at the position of the 
perturber x — y = 0, z — vt, 

d x (p F = 8 X ^ F = d y (j) F — <3j,V F = EnotDF , (59) 

as expected from the cylindrical symmetry of the problem, while 
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the t and z gradients, evaluated at the position of the perturber x 

y = 0, z — vt, are 



9 Z V F = d z <f> f 



_dt£ 

V 

0<*- 



+ 1/ 2 + (z'-,3) 2 ] 3 / 2 



X (1 + £ Yukawa) + £not DF ■ 



(60) 



[Note that this expression for dtip F is obtained by taking the deriva- 
tive of equation J57| > with respect to t, transforming the derivative 
with respect to t acting on 8n/n into a derivative with respect to z' 
using equation d43l >. integrating by parts and finally transforming 
the derivative with respect to z' into one with respect to z.] The 
integral in equa tion d60t can be evaluated using equation d43 1 as in 
Ostriker ( 1999), and is 



<9zV> = d z (/> 



d t <f> ¥ 



(61) 



4tt(p + p)M-y(l + v 2 



I x [1 + 0(C/\j)} + e no t 



DF , 



where 



I = 



iln 



l+M 
1-M 



-M 



ifM < 1. 



iln(l-^)+ln(^-) \fM>l 



(62) 



and we have made the assumptions that \c s — v\t exceeds the cutoff 
r m i n and that \c 3 + v\t is smaller than C. 

Inserting equation l |59t into equation 1471 . one immediately 

finds 



(Spert) F — («pert) F — Enot DF 

while using equation J61t in equation (47) gives 

/-« \ _ 47v(p + p)Mr\l + v 2 ) 2 . 



(63) 



'■port J F 



(64) 



(5pert) F 



x [1 + 0(C/\j) + £>(M/r min )] + e not df 
4tt( P + p)Mj 3 (l + v 2 ) 2 

x [1 + 0{C/\j)+0(M/r min )] + e notDF 



(65) 



where / is defined by equation l |62t . Note that we have restored the 
relative error C(ei, £2) due to the linearization of the equations of 
the previous section: this gives rise to the error 0(M/r m i n ) ap- 
pearing in equations l |64t and i65\ . 

Performing a boost we can calculate the change of 3- 
momentum in the rest frame of the perturber due to the gravita- 
ti onal interaction with the fluid, so as to compare with the results 
of |Petrichet"ai] ( ll989h : 



ert 



dr 



M"f [(5p Crt ) F - v (apert)i 



(66) 



4n(p + p)M 2 -y 2 (l + v 2 



■Ix[l + 0(C/Xj) + C(A//r min )] 



+ £not DF , 

■fx) > 
pert 



J -(a) 



dv (v) 



— Enot DF • 



(67) 



dr I \ dr 

/ F \ / F 

Note that the relative errors 0(C/Xj) and C(M/r m i n ) are negligi- 
ble - the former because the fluid is not self-gravitating and the lat- 
ter because the effective cutoff radius r m i n is large compared with 
M - whereas Snot df in general is not negligible. However, e no t df 
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Figure 1. The relativistic correction factors 7 2 [1 + (v/c) 2 ] 2 , multiplying 
the Newtonian drag for straight-line motion and the tangential Newtonian 
drag for circular motion, and 7 3 [1 + (v/c) 2 ] 2 , multiplying the radial New- 
tonian drag for circular motion, are plotted as functions of the velocity v 
of the perturber relative to the fluid. Note that velocities v ~ 0.8 can be 
obtained for a perturber orbiting around an accreting SMBH in the opposite 
direction with respect to the accretion flow ( Barausse 2007). 



represents the standard force acting on the perturber because of the 
gravitational interaction with the fluid, and it can be computed sep- 
arately if the global structure of the system is known. In particular, 
£not df = if the medium is distributed in a spherically symmetric 
fashion around the perturber. 

The relativistic correction factor 7 2 (1 + v 2 ) 2 appearing in 



equation d66| > is plotted as a function of the velocity v in Fig. 
[T] Note that, for A4 3> 1 and vt —> r max , equation j66\ , and 
in particular the co rrection factor, agrees with equation (B45) of 
|Petrichetai] ( ll989h . 



4 CIRCULAR MOTION 

Let us now consider the case of a perturber moving on a circular 
orbit of radius R with velocity v — flR. Such an orbit is clearly 
not allowed in a Minkowski background, unless there is an ex- 
ternal non- gravitational force keeping the perturber on a circular 
trajectory. In astrophysical scenarios we are interested instead in 
a perturber maintained in circular motion by gravitational forces. 
In this case, the background spacetime is necessarily curved: one 
can think of a circular orbit around a Schwarzschild black hole 
with mass A/bh surrounded by a tenuous fluid at rest. However, 
if the perturber is sufficiently far from the central black hole {i.e., 
if R S> Mbh) one can approximately consider the metric as given 
by equation d 1 5b (i.e., Minkowski plus the perturbations produced 
by the presence of the fluid and of the perturber) and neglect the 
corrections 0(Mbu/ R) due to the presence of the central black 
hole. This treatment is clearly not completely satisfactory, because 
orbital velocities become relativistic only close to the central black 
hole [in fact, v ~ (Mbh/ R) 1 ^ 2 ], but we argue that it may not be 
such a bad approximation as it might seem. 

Indeed, if one uses Fer mi normal coordinates comov i ng wi th 
the perturber [see for instance lMisner. Thome. & Wheeled dl 9731) 1 ■ 
all along the trajectory the metric can be written as Minkowski 
plus perturbations produced by the fluid and the perturber, the 
curvature of the background introducing just corrections of order 
0(r /Mbh) 2 (r being the spatial distance from the perturber). Be- 
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cause the wake can extend o ut to distances of order R from the 
perturber dKim & Kimll2007h . it will eventually feel the curvature 
of the background unless R <C A/bh. However, the part of the 
wake giving the largest gravitational attraction to the perturber will 
be the closest to it, and this part will experience an approximately 
flat spacetime. Reasoning in the same way, we can argue that our 
treatment should be approximately applicable also to a perturber 
moving on a circular orbit in a fluid which is moving circularly in 
the same plane as the perturber (e.g., a perturber moving inside an 
accretion disc), provided that the velocity v = QR of the perturber 
is taken to be the velocity relative to the fluid. 

Considering therefore a Minkowski background spacetime, 
one can proceed as in the previous section, and equation J34t be- 



gradient, instead, we have 



(72) 



(p + p) J d*x 

X (1 + £ Yukawa) + £ not DF 

4ir(p + p)My(l + v 2 )R 



,%(x',t = 6/S}) (x-x')-dex 



Is X [l+e>(£/A./)]+£notDF , 



where Ig is given by 



M 2 
4% 



]3 "/ 

a x 



V{x',t = 9/n) r'sm(0' -0) 
[1 + z' 2 + f' 2 - 2f' cos(6> - 6")] 3/2 



(73) 



(fl? _ c , y2) ^ = 4nM 7 (l + v 2 ) s{r R)s{z)5{e 

+ 4w[p(l + 20) + 3p(l - 2-0)] + 4tt(1 + 3c 2 s )(p + p) — , 

(68) 



where we have introduced a system of cylindrical coordinates 
(r, 9, z) such that the motion of the perturber takes place at z = 0, 
r — R. The solution to this equation is ra t her co mplex, but has for- 
tunately been worked out bv lKim&Kin]j2007h . For our purposes, 
proceeding as in the previous section we can simply write it as 



-V(x,t)x(l + £3)+0{C/\.,y , (69) 



n Ret 



where ej ~ 0(C/\j) and t he weight-function T>(x,t), whose 
detailed form can be found in iKim & Kiml J2007h . defines the re- 
gion of influence which sound waves sent off by the perturber do 
not have time to leave. From the plane-symmetry of the problem, 
it is clear that T)(x, y, z, i) = T>(x, y, — z, i) and ej(x, y, z, t) = 
£j (x, y, —z,t). Moreover, from the gradient of equation i6S\ it also 
follows that 



d t V = -fi d V 



(70) 



If we are again concerned with the force exerted by the fluid, which 
includes dynamical friction effects, rather than with the accretion 
drag or the self-force 0, we can restrict our attention to the met- 
ric perturbations 4> F = ip ¥ generated by the fluid, which are again 
given by equation l |49t . Using again the Green's function of the 
Yukawa operator and evaluating at the position of the perturber 
(r — R, 6 — ttt, z — 0) we easily get 



8 z <j) ¥ = d z 1p F = Snot DF 



(71) 



(from the plane symmetry of the function T>). For the azimuthal 



B Note that, differently from the case of straight-line motion, even the dis- 
sipative part of the self-force is now non-zero, as can be seen from the fact 
that the perturber loses energy and angular momentum through gravitational 
waves (cf. the quadrupole formula). S elf-force calcu lations, as already men- 
tioned, require different techniques iPoisson 2004) and can be performed 
separately. 



(a hat denotes quantities scaled by the radius of the orbit: x' = 
x' /R, r' — r' I R, z = z'/R). Similarly, for the radial gradient 
we obtain 



9 r -V ,F = d r (p F = 



(74) 



, %(x',t = 0/n) (x-X^drX 

(P + P) / d x r- — — t, 



X (1 + EYukawa) + Enot DF 

47t(p + p) Af7(l + v 2 ) 



I r x [l + 0(C/Xj)]+e not T>F 



where is given by 



Ir = 



M 2 

47T 



z ^ V{x',t = 9/Q) [r' cos(0-9') - 1] 
X [1 + z' 2 + r' 2 - If' cos(0 - 0')13/2 



(75) 

N ote that the integral s Ie and I r have been calculated numerically 
in lKim&Kiml (2007). They are functions of the coordinate 9 of the 
perturber, which is thought to vary in an unbound range to count 
the number of revolutions, or equivalently they can be thought of 
as functions of time (t — 9/Q). Fortunately, though, steady state 
values for these integrals are reached in times comparable to the 
sound crossing-time R/c a or within one orbital period: fits to the 
numerical results for these steady state values i n the case in which 
R > r min and C > (20 - 100) R are given by dKim & Kiml2007h 

< M 2 10 3 B1M - 4 22 , for M < 1.1, 



0.5 In [9.33M 2 {M 2 -0.95)] 
for 1.1 < M < 4.4, 

0.3 M 2 , forX > 4.4, 



(76) 



and 



0.7706 In, l 0004 _ u:) | N -vi 
for M < 1.0 , 



1.4703X, 



(77) 



ln[330(i?/r mi „)(X - 0.71) 5 ' 72 X" 9 ' 58 ] 
for 1.0 < M< 4.4, 

ln[(i?/r min )/(0.11M + 1.65)], 
for M>4A. 



These fits are accurate within 4% for A4 < 4.4 and within 16% for 
M > 4.4. 

Using equations l |7U , l |72t and l |74| l in equation l |47t and trans- 
forming to cylindrical coordinates, for the acceleration produced by 
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the gravitational interaction with the fluid we easily get 

( -t \ Mp + p)Mj 3 (i + v 2 ) 2 

( a pert] F ie 



(78) 



f-e \ 

I a pert I 



x [1 + 0{C/Xj) + 0{M/ 7"min)J ~\~ £not DF 

4 7 r(p + / 9)M 7 3 (l+w 2 ) 2 T 

= le 

F Hv 

x [1 + 0{C/Xj) + 0(M/ ) + Snot DF 



(79) 



^pert )p — 



4iv(p + p)My 3 (l + v 2 ) 2 



It 



x [1 + 0(C/Xj) + 0(M/r min )] + e not , 



(80) 



(81) 



[The error 0(M/r m i n ) comes about because the equations that we 
have solved are linearized and are therefore subject to an "intrinsic" 
error 0(ei, £2).] 

Finally, in order to compute the change of 3-momentum due 
to the gravitational interaction with the fluid in the rest frame of 
the perturber, it is sufficient to project the 4-force M(dp elt )F onto 
a tetrad comoving with the perturber, i.e., e( t ) = u^ ert d/dx^ = 
7(0/9* + Sid/ 09), e w = j(vd/dt + l/rd/80), e (r) = d/dr 
and er z -\ = d/dz. Using equations J78b — J8 1 b one then easily gets 



i^pcrt 

dr 



47r(p + p)M 2 7 2 (l +v 2 ) 2 



x [1 + 0(C/Xj) + C(M/r min )] + e not DF , (82) 

dp ( p r L \ _ MP + p)M 2 J 3 (l+v 2 ) 2 T 
dr 



x [1 + 0{C/Xj) + 0(M/ ) + 

£not I 



dv (z) 

"-Fpcrt 

dr 



— Enot DF • 



(83) 
(84) 



As in the case of straight-line motion, the relative errors 0(jC/Xj) 
and 0(M/r m i n ) are negligible, because the fluid is not self- 
gravitating and because the effective cutoff radius r m i n is large 
compared with M, whereas £ no tDF in general is not negligible, 
although it is exactly zero if the medium is spherically symmetric 
around the perturber. The relativistic correction factors 7 2 (l + f 2 ) 2 
and 7 3 (1 + « 2 ) 2 appearing in equations d82t and ( 1831 ) are plotted 
as functions of the velocity v in Fig.Q] 



5 CONCLUSIONS 

We have studied the drag experienced by a massive body because 
of the gravitational interaction with its own gravitationally-induced 
wake, when it is moving along a straight-line or a circular orbit 
at relativistic speed v relative to a non self-gravitating collisional 
fluid in a flat or weakly curved background spacetime. Thanks to 
a suita ble choice of ga uge, w e could explo i t the Newtonian anal- 
ysis of lOstrikej \ 19991) and of iKim & Kim ( 2007) to simplify our 
calculations. We find that their results remain valid also in the rela- 
tivistic case, provided that the rest-mass density is replaced by p+p 
(p and p being the pressure and energy density of the fluid) and a 
relativistic multiplicative factor is included. This factor turns out 
to be 7 2 [1 + (v/c) 2 ] 2 in the straight-line motion case and for the 
tangential component of the drag in the circular motion case, and 
7 3 [1 + (v/c) 2 ] 2 for the radial component of the drag in the circular 
motion case. 

Although our analysis strictly applies only to a fluid in a flat 



spacetime (in the case of straight-line motion) or a weakly curved 
one (in the case of circular motion), we have argued that our re- 
sults are suitable at least for a preliminary study of the effects of 
an accretion disc on EMRIs. Although our results are not expected 
to change the standard conclusion that the gas accreting onto the 
central SMBH does not si gnificantly affe ct EMRIs in the case of 
"normal" Galactic Nuclei dNaravanll2000h , they could play a role, 
under certain circumstances, in the case of higher density environ- 
ments like Active Galactic Nuclei (quasars, Seyfert Galaxies, etc.). 
An investigation of this scenario, in which the accretion is modelled 
by a t hick torus, will be presented in a subsequent paper dBaraussel 
l2007h . 
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